"Locally homogeneous turbulence" Is it an inconsistent framework? 
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In his first 1941 paper Kolmogorov assumed that the velocity has increments which are homoge- 
neous and independent of the velocity at a suitable reference point. This assumption of local homo- 
geneity is consistent with the nonlinear dynamics only in an asymptotic sense when the reference 
point is far away. The inconsistency is illustrated numerically using the Burgers equation. Kol- 
mogorov's derivation of the four-fifths law for the third-order structure function and its anisotropic 
generalization are actually valid only for homogeneous turbulence, but a local version due to Duchon 
and Robert still holds. A Kolomogorov-Landau approach is proposed to handle the effect of fluctu- 
ations in the large-scale velocity on small-scale statistical properties; it is is only a mild extension 
of the 1941 theory and does not incorporate intermittency effects. 



The concept of homogeneity, which goes back to 
Kelvin, 1 applies to flows whose statistical properties are 
invariant under (space) translations. For translation- 
invariant equations, such as the Burgers or Navier-Stokes 
equations and in the absence of boundaries it is obvious 
that homogeneity is dynamically preserved. More pre- 
cisely, with a homogeneous initial velocity field uq, as 
long as the solution is unique, it will also be homoge- 
neous. 

A less restrictive condition is that of incremental ho- 
mogeneity: we now demand only that the statistical 
properties of the increments u(x + r) — u{x) be invari- 
ant under translations, that is depend only on r but 
not on x. In the time domain for which homogeneity 
is called stationarity, a well-known example of an incre- 
mentally stationary random function which is not sta- 
tionary is the Brownian motion curve. As pointed out 
by Shiryaev^ Kolmogorov published in 1940, one year 
before his famous publications on turbulence, two prob- 
ability papers. One discusses among other things incre- 
mental homogeneity 3 and the other one "Wiener spirals" , 
self-similar Gaussian random functions now called frac- 
tional Brownian motions— He was dealing with random 
functions whose energy spectrum E{k) has an infrared 
divergence, that is J °° E{k)dk diverges at k — 0. For ex- 
ample the Brownian motion has E(k) cx k~ 2 . Turbulence 
in fluids and plasmas is replete with infrared divergent 
power-law spectra, for example, the Kolmogorov fc~ 5 / 3 
spectrum. In such cases the variance of the velocity in- 
crement is finite but the variance of the velocity itself 
is infinite. Of course, upon introduction of an infrared 
cutoff fco > 0, homogeneity is recovered. For example 
the homogeneous Ornstein-Uhlenbeck Gaussian process 
which has the spectrum E(k) = C '/(fco + goes over 
into Brownian motion as fco — * 0. 

In his first 1941 paper— (K41) Kolmogorov pointed 
out the need to work with the increments of the ve- 
locity. However, rather than assuming that the tur- 
bulence is incrementally homogeneous, he assumed lo- 
cal homogeneity. Specifically he took a space-time ref- 
erence point p(°) = (t*(°), i(°)) and, for an arbitrary 



space-time point P, considered the space-time increments 
w(P) = u(P) — tto where Uo — u(P^) and then defined 
local homogeneity as the property that for any choice of 
n points Pi, P2, P„, the joint distribution of the 
n increments should be independent of both p(°) and 
of UQm^ In other words, Kolmogorov was assuming not 
only incremental homogeneity and incremental station- 
arity, he was also assuming independence of the velocity 
at the reference point. 

As stressed by Hillj£i£ various definitions of local ho- 
mogeneity have been used in the literature. For example 
Monin and Yaglom& mostly take it to be equivalent to 
incremental homogeneity but, of course, also quote Kol- 
mogorov's full definition. It is not clear why Kolmogorov 
felt it necessary to introduce the additional assumption 
of independence of the velocity at the reference point. 
It is conceivable that when he wrote his first 1941 tur- 
bulence paper he already knew that he would need this 
assumption to avoid extra terms appearing otherwise in 
the derivation of the four-fifths law— In this paper we 
shall discuss successively the problems raised by incre- 
mental homogeneity and those raised by the assumption 
of independence of u$. 

First we observe that it is far from obvious that incre- 
mental homogeneity is consistent with the quadratic non- 
linear dynamics of the Burgers and Navier-Stokes equa- 
tion. Indeed the square of an incrementally homogeneous 
function has no reason to share this property because the 
increments of the square are not the squares of the incre- 
ments. 

It is easy to show the dynamical inconsistency of incre- 
mental homogeneity using the one-dimensional Burgers 
equation 



d t u + d x — = vd xx u. 



(1) 



We shall work with the solution in the limit of vanishing 
viscosity which has a simple representation: writing u = 
—d x ip, we have 



ip{x,t) 
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where *0(x,O) is the initial (Lagrangian) velocity poten- 
tial (see, e.g. Ref.[l(|. Eq. © can be implemented nu- 
merically in a very efficient way using the Fast Legendre 
transform algorithm^! A particularly interesting initial 
condition which has been much studied (for example be- 
cause it arises in cosmologj*ii*iSiii) is to take for u(a, 0) 
the bilateral Brownian motion curve passing through the 
origin (see Fig.^l. For a > it is defined as a Gaussian 
random function with zero mean and with correlation 
function (u(a, 0)u(a' , 0)) = inf(a, a'). For a < it is 
defined as u{— a,0) where w(a,0) is another realization 
of u(a, 0) independent of u(a,0). With such choice of a 
self-similar initial condition it is easily shown that the 
statistical properties of the solutions at any two given 
non-vanishing times are related by a simple change of 
scaled Thus no generality is lost by assuming t = 1, as 
we shall do here. 

In practice u{a, 0) is generated as a bilateral ran- 
dom walk (with steps ±V Aa ) on a regular grid with 
mesh Aa = 10 -3 , covering the Lagrangian interval 
[— Liag, +Li ag ] where L\ ag = 1 3 . The solution is calcu- 
lated in a much smaller Eulerian interval [— L cu i, +L cu i] 
with L ou i = Liag/100 = 10 to ensure that the probabil- 
ity to have a Lagrangian antecedent outside the interval 
[— Liag, +iiag] is negligibly small. It is important not to 
take initial conditions just from a Fourier series since this 
would make the initial condition periodic, thereby loosing 
incremental homogeneity in favor of homogeneity, which 
is dynamically preserved. Periodicity, assumed in pre- 
vious numerical studies^i thus cannot reveal what we 
see here. Fig. |21 shows the second-order structure func- 
tion ((u(x + r) — u(x)) ) obtained after averaging over 
10 5 realizations of the initial conditions. It is here plot- 
ted against the absolute position x at fixed separation r. 
A strong dependence on x reveals that incremental ho- 
mogeneity does not hold, except at very large x where the 
structure function becomes translation- invariant. The 
validity of the numerical results was tested by increasing 
by one order of magnitude the length of the Lagrangian 
interval and decreasing the number of realizations so that 
the computation remains manageable. Exactly the same 
picture as before is obtained, but of course with noisier 
statistics. 

A proof of the loss of incremental homogeneity in the 
temporal evolution from an initially incrementally homo- 
geneous velocity can be given for the case of smooth ini- 
tial data. The general idea is that as long as there are no 
shocks the solution can be Taylor expanded in time. The 
second-order structure function can then be calculated 
perturbatively and shown not to be translation-invariant. 
Actually for Gaussian initial data some of the realizations 
will develop shocks after very short times but their con- 
tribution to the structure function can be bounded by 
terms which are exponentially small in 1/t 2 . Because of 
its perturbative nature it is likely that the proof can be 
extended to the case of the 3-D Euler or Navier-Stokes 
equations. The proof is quite long and technical; its de- 
tails do not belong here. 




FIG. 1: A snapshot of a Brownian initial velocity (thin line) 
and the corresponding solution to the Burgers equation at 
time one (thick line); notice the proliferation of shocks. 




FIG. 2: The second-order structure function for various sep- 
arations r, plotted vs the absolute position x. Notice the 
strong x-dependence, evidence that incremental homogeneity 
does not hold. At large x the graphs become flat (not shown). 



Since incrementally homogeneous random functions 
may be obtained by infrared limits of homogeneous ran- 
dom functions for which homogeneity persists under 
Burgers or Navier-Stokes dynamics, it is of interest to un- 
derstand how homogeneity and incremental homogeneity 
are lost when the limit is taken. For concreteness we 
shall consider the Burgers equation with an initial veloc- 
ity which is the Gaussian Ornstein-Uhlenbeck process 
(OUP). The argument given below is actually very gen- 
eral and can be applied mutatis mutandis to the Navier- 
Stokes equation. The OUP u(x, 0; fco) has the spectrum 
1/(^0 + k 2 ). Hence its correlation length L ~ 1/fco and 
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its variance (u 2 (x, 0; fco)) ~ V^o §° t° infinity as fco — > 0. 
The Brownian motion process must be equal to zero at 
the origin. It is obtained from the OUP by first sub- 
tracting the random velocity u(0, 0; fco) and then letting 
fco — > 0. We are here emphasizing "random" because 
there is a similar stress put by Kolmogorov on the first 
page of his first 1941 paper. It is clear that the effect 
of the initial subtraction of a uniform velocity will be a 
(random) Galilean transformation. Hence, denoting by 
u(x,t;ko) the solution of the Burgers equation with the 
initial condition M(x,0;fcg), the solution of the problem 
with u(0, 0; fco) subtracted is 

v(x, t; fco) = u(x + tu(0, 0; fco), t; fco) — u(0, 0; fco). (3) 

As fco — > the function v(%, t;ko) tends to the solution 
with Brownian initial condition which was determined 
numerically above. Inspecting the r.h.s. of © we note 
that if we had only the subtraction of w(0,0;fco) the so- 
lution and its limit would remain incrementally homoge- 
neous. However, we also have a dependence on u(0, 0; fco) 
inside the spatial argument. Under the assumption of lo- 
cal homogeneity, the r.h.s. of J3J would be independent 
of the velocity u(0, 0; fco) at the reference point (here the 
origin) and thus v would be incrementally homogeneous. 
It is this violation of local homogeneity (after subtrac- 
tion of u(0,0;fco)) which spoils both homogeneity and 
incremental homogeneity. Local and incremental homo- 
geneity can only hold in an asymptotic sense at points 
sufficiently remote in space-time from the reference point 
to be only weakly influenced by its velocity. 

At the very beginning of his third 1941 turbulence 
paper on the four-fifths law Kolmogorov 9 assumed lo- 
cal homogeneity^ 5 How is the four-fifths law affected 
by what we just discussed? Monin and Yaglom^ real- 
ized that there may be problems in deriving the four- 
fifths law (and its anisotropic generalization) using lo- 
cal homogeneity rather than full homogeneity. Hence 
they proposed first (on p. 395) to consider the locally 
homogencous/isotropic turbulence as being embedded in 
a large-scale homogeneous turbulence. Then (on pp. 401- 
403) they derived the same results for locally homo- 
geneous/isotropic turbulence using a quasi-Lagrangian 
method with a frame attached to a particular fluid par- 
ticle. Actually, as we have seen, neither incremental 
nor local homogeneity can hold for turbulence which is 
initially not strictly homogeneous, except in an asymp- 
totic sense which would not help in deriving a four-fifths 
law. Even in a putative incremental homogeneous turbu- 
lence, attempts to derive a four-fifths law lead to an extra 
non- vanishing term as shown in Ref. 14j (second term in 
Eq. (39)) and in Ref . (second term in Eq. (21)). This 
term is associated to the advection of products of two 
velocity differences by the half sum. Of course, the extra 
term vanishes if homogeneity (or local homogeneity) were 
to hold. We cannot rule out that some conditional form 
of homogeneity holds which allows a dependence on the 
reference velocity but is consistent with the four-fifths 
law. 



When the turbulence is initially not homogeneous we 
cannot write a standard four-fifths law but there is still a 
relation between the local fluctuating energy dissipation 
s(x) and the velocity increments, due to Duchon and 
Robert In it simplest form it reads 

e(x) = lim ^ f dfl \u(x+r)— u(x)\ 2 (u(x+r)— u(x))-r, 

(4) 

where d£l is the element of solid angle spanned by the 
unit vector r/r. This relation, which holds only after 
the limit v — ► has been taken, is between non-averaged 
quantities. In unpublished notes by Onsager (cited by 
Eyink and Sreenivasani&) the same relation appears in 
averaged form^& The Duchon-Robert formula is impor- 
tant because it is an exact consequence of the dynam- 
ical equations and also because it can be viewed as an 
infinitesimal version of Kolmogorov's 1962 Refined Sim- 
ilarity Hypothesis. 17 

We have seen that Kolomogorov's assumption of in- 
dependence of statistical properties of increments on the 
velocity at the reference point holds at best in an asymp- 
totic sense, far from the reference point. More generally, 
one can ask if the small-scale properties of turbulent flow 
near a given point depend much on the instantaneous 
velocity uo at this point. There is good evidence that it 
does. For example, there are recent experimental^ and 
then numericaiiS data showing a strong dependence of 
the variance of the fluid particle acceleration, conditioned 

9 /2 

on wo, which is found to vary roughly as u . Indeed as 
shown in Ref. |2(], large velocity and (extremely) large ac- 
celeration events occur in a correlated way near vortex 
filaments. 

There is nothing surprising to find dependence on uo 
in homogeneous and isotropic turbulence: if in such a 
flow we observe a high fluctuation of the local velocity 
we also expect a high fluctuation in the local Reynolds 
number, a sharp decrease in the dissipation scale, etc. 
Some of the dependence on uq can be captured by adapt- 
ing an old argument of Landau (see Ref. l2ll Sec. 6.4.2). 
In its original form, Landau's argument assumed that 
the mean dissipation e which appears in the Kolmogorov 
1941 expression for the structure functions, can actually 
be taken as a quantity presenting large-scale fluctuations. 
As pointed out in Ref. |52 arguments a la Landau have 
the potential of interlinking large-scale and small-scale 
phenomena. For this, one just uses the Kolmogorov 1941 
predictions in terms of the large scale velocity uq, the in- 
tegral scale L and the viscosity (when needed) and then 
one allows fluctuations in uq*L For example, as pointed 
out in Ref. l23l from the Heisenberg-Yaglom theory which 
predicts that the typical acceleration of fluid particles is 
proportional to Ug^ 4 , one can then infer that its condi- 
tional variance goes as Uq^ 2 . If Mg is assumed to have a 
Gaussian distribution, the PDF of the acceleration will 
involve a stretched exponential. 1 - Such non-Gaussianity 
should not be mistaken for intermittency. It is just a con- 
sequence of a slight extension of the Kolmogorov theory, 
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which might be called Kolmogorov-Landau theory. 

A delicate issue concerns the effect of it on small-scale 
velocity increments and in particular on the longitudinal 
and transverse second-order structure functions. This 
question is now being revisited by various colleagues. 
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